CHAPTER IX 


ADAPTIVE ATTITUDE CONTROL OF THE 
CREW LAUNCH VEHICLE 

9. 1 Introduction 

Classical linear control laws have been developed for the NASA’s Crew Launch Vehicle 
(CLV). Preliminary design of the Ares-I flight control system has shown that classical 
control theory is sufficient to meet the stability and performance requirements of 
the CLV[10]. Due to the uncertain nature of a launch vehicle, imprecise knowledge of 
system parameters and the flexibility of the CLV rocket during its ascent phase creates 
concern that unforeseen instabilities may develop. Time varying factors such as fuel 
consumption and mass reduction, combined with a wide range of aerodynamic and 
dynamic interactions including payloads, propulsion, inertia, and dynamic pressure, 
can affect the overall dynamics of the vehicle during its flight. 

Studies over the past few years have suggested that adaptive control techniques 
can potentially be beneficial to the flight control system in terms of robustness and 
safety [122], Adaptive control may not only provide a higher level of nominal perfor- 
mance, but can also accommodate a greater degree of uncertainty, including active 
vibration suppression of uncertain flexible modes[129] and accommodation of par- 
tial failures in the flight control system. More specifically, it has been shown that 
adaptive control can increase the robust performance of the CLV[88]. However, the 
nature of an adaptive control signal raises concern for implementation on a manned 
launch vehicle. Since model reference adaptive control laws are inherently high band- 
width, high frequency control effort maybe generated that can destabilize unmodeled 
dynamics and exceed actuator capabilities. Moreover, if the uncertainty does not 
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satisfy matching conditions, the adaptive control law can degrade tracking and even 
destabilize the plant. 

This major disadvantage of adaptive control stems from the fact that it lacks 
an accepted means of quantifying the behavior of the control signal a priori. These 
measures are required in order to certify flight control systems of piloted and pas- 
senger bearing vehicles. Hence, most adaptive control laws require a more extensive 
verification and validation process due to the time varying and nonlinear manner in 
which its gains are adapted since the control effort could be beyond the limits of the 
system. From this prospective, it is highly desirable to limit the frequency content 
of an adaptive control signal. Classic and robust control offer natural frameworks 
for achieving frequency limited signals. The H ^- NMA Architecture [90] allows one 
to also achieve frequency limited control signals by combining aspects from robust 
control theory and adaptive control theory into a single framework. This makes it 
possible to limit the frequency content of the adaptive control signal in an algorithmic 
manner. 

In this chapter, attitude control of the CLV is accomplished using two decou- 
pled Hoo-NMA state feedback architectures designed to maintain the design level of 
tracking performance in the presence of disturbances and parametric uncertainties. 
Emphasis is placed on a minimal order adaptive law in order to see if a low band- 
width, low gain, and reduced order state feedback control law can offer performance 
improvement. This new control architecture merges ideas from LAo control theory 
and adaptive control theory to achieve band limited control signals. This represents 
a different approach than previous approaches based on an high order output feed- 
back a-modihcation adaptive law[88]. To show the viability of the method, a high 
fidelity simulation of the CLV(called SAVANT) is used in conjunction with the actual 
decentralized classical control design used on the CLV to compare the nominal perfor- 
mance against the performance of the CLV with an augmented adaptive law. Rigid 
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body, aerodynamics, gravity, sloshing, engine inertia effects, mass change, actuator, 
and clastic body models(among other things) are included in the simulation. The 
nominal controller consists of three independently designed controllers for yaw, pitch 
and roll attitude. The effect of structural modes are compensated for with using gain 
and phase stabilization filters in the pitch and yaw channels. The roll control law is 
a nonlinear bang-zero-bang design. The presence of the bang-zero-bang control law 
requires a special modification to the adaptive law called control hedging in order to 
ensure that the nonlinear nature of the roll control law if properly accounted. This 
has been shown to work well [88]. However, the roll control channel in the work is 
ignored in the adaptive design due to design restrictions from NASA. In implemen- 
tation, the original control design is not modified. The adaptive control law simply 
augments the nominal control law. This facilitates switching the adaptive control law 
on in case of degraded nominal control performance. 

Results examine the degree to which the nominal control design can be improved 
by adding an adaptive element. To measure the degree of improvement, the Worst- on- 
Worst(WoW) Monte-Carlo dispersion cases are compared. The WoW cases capture 
the combination of the worst possible uncertainties (i.e. dispersions) occurring simul- 
taneously. Simulation results show that the adaptive control law always improves the 
performance of all of the stable WoW cases. For the WoW cases that are unstable 
with only the linear control law, the adaptive control algorithm is able to maintain 
acceptable tracking performance in almost all of the cases. 

9.2 Vehicle Model and Dynamics 

The Ares I CLV is a two-staged, serially connected rocket with the Orion crew ex- 
ploration vehicle located at the top. The launch vehicle’s first stage consists of a 
single, five-segment reusable solid rocket booster, and the second or upper stage is 
propelled by a main engine fueled with liquid oxygen and liquid hydrogen. The vehicle 
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configuration is shown in figure 9.1. 


♦ ~25-mT payload capacity 

♦ 2-Mlb gross liftoff weight 

♦ 309 ft in length 


Upper Stage 

• 280-klb Liquid Oxygen/Liquid Hydrogen 
(LOX/LH 2) Stage 

• 5.5-m Diameter 

• Aluminum-Lithium (Al-Li) Structures 

• Instrument Unit and Interstage 

• RCS / Roll Control for First Stage flight 

• CLV Avionics System 


First Stage 

Derived from Current Shuttle Reusable 
Solid Rocket Motor/Booster (RSRM/B) 
Five Segments/Polybutadiene 
Acrylonitride (PBAN) Propellant 
Recoverable 
New Forward Adapter 



Upper Stage Engine 

• Saturn J-2 Derived Engine (J-2X) 

• Expendable 


Figure 9.1: Ares I Crew Launch Vehicle 


Ares I has two vital missions; lifting astronauts up to the International Space Sta- 
tion and achieving an in orbit rendezvous with the Ares V Earth departure stage at 
low Earth orbit for a mission to the moon. During the first two and a half minutes of 
flight, the first stage booster powers the vehicle to an altitude of about 38 miles and 
a speed of Mach 5.9. After its propellant is consumed, the solid rocket booster sepa- 
rates. The upper stage engine is then ignited and powers the Orion spacecraft. After 
reaching an altitude of 83 miles, the upper stage separates and the Orion spacecraft 
completes its trip to a circular orbit of 185 miles above the Earth using its service 
module propulsion system. 

9.2.1 CLV Model Description 

The CLV model employed in this study is called Savant and was developed in a joint 
effort between bD Systems and NASA Marshall. The model is described in Betts [9] 
and it contains simulated rigid body, aerodynamics, gravity, sloshing, engine inertia 
effects, mass change, actuator, and elastic body models. Many of these effects can be 
turned on or off in the simulation environment. 
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Since the Ares I CLV possesses the characteristics of a long and slender body, its 
flexibility should be considered in the control law design. In the structural modeling 
part, modal frequencies, displacement, and rotation are given from a Nastran (FEM 
solver) solution and is used to model the interaction effects between the vehicle flex- 
ibility and the other dynamic models. Lateral vibration is the dominate vibration 
mode. It is important to consider the effect of this vibration in control system design 
due to the modal frequencies being near the control bandwidth. The vehicle’s elastic 
motion can be conveniently expressed in terms of frequencies and mode shapes of a 
free-free beam structure. Because of the axial symmetry of Ares I launch vehicle, 
two identical modes exist in the lateral bending. Table 9.1 gives a summary of the 
approximate dominate bending mode frequencies at launch. The actual modal fre- 
quencies are time varying and change significantly through out ascent. Figure 9.2 
shows the dominate vehicle mode shapes. 


Table 9.1: Bending Frequency Table 


1st Bending Frequency 

2nd Bending Frequency 

3rd Bending Frequency 

6.0 [rad/sec 

14.2 [rad/sec] 

27.2 [rad/sec] 


The aerodynamic model contains three parts; the environment model, the aerody- 
namic coefficients, and the force and moment generation algorithm. The time history 
of the mass properties including the total vehicle weight and the center of gravity 
location is stored in the simulation model in the form of look-up tables. Multiple 
gimbal actuator models are available for simulating actuator limits including a 3rd 
order model or a high fidelity simplex model. The fuel sloshing model simulates the 
effect of fuel sloshing as point masses connected to the rocket body as shown in figure 
9.3 to the LOX and LH2 upper stage tanks by a spring and damper in the lateral (Y, 
Z) direction. The vertical position of slosh point mass is function of the liquid level 
in the tank. Vehicle separation is modeled as an instantaneous loss of mass. The 
atmosphere is simulated using the US76 atmosphere model and the gravity model 
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0.02 



(a) First Bending Mode 



(b) Second Bending Mode 



(c) Third Bending Mode 

Figure 9.2: Visualization of the first three structural mode shapes 
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Figure 9.3: Schematic of elastic vehicle with sloshing point mass and engine inertia 

used includes Earth oblateness effects without considering abnormality or vertical 
deflection data. 

9.2.2 Nominal Control Law 

The nominal control law tracks quaternion guidance commands. At each time instant, 
the quaternion command and the vehicle’s attitude quaternion are used to generate 
an attitude error signal. These errors are suppressed using two gain scheduled PID 
control laws for the pitch and yaw degrees of freedom and a phase plane roll controller 
for the roll degree of freedom. One restriction of the roll RCS is that the actuator 
only fires when the roll error of the CLV is greater than a defined parameter (a bang- 
zero-bang control law). Each axis is assumed independent so each control channel 
is designed separately. The effects of structural modes are gain and phase stabilized 
using a combination of gain scheduled low pass and notch filters. Gains are scheduled 
based on the altitude. The nominal control architecture is shown in figure 9.4. 

The error signal to the nominal control law is generated based on a quaternion 
command signal. A comparison between the command and the current vehicle quater- 
nion state generates an error angle state that is used for feedback. This error angle 



Figure 9.4: Diagram of nominal control system 


state is the incremental Euler angle rotation from the current system state to the cur- 
rent command. To define this error angle signal, suppose that there are two frames, 
frame 1 and frame 2. One can describe the rotation from frame 1 to frame 2 in terms 
of an axis and an angle to rotate around that axis. This representation can be used to 
define the quaternions used in this chapter. The following definition of the quaternion 
vector based on an axis-angle formulation is used. 


Qaxis— angle 


cos(£/2) 
a x sin(£/ 2) 


(9.1) 


a y sin(£/ 2) 
a z sin(£/ 2) 

where a is a unit vector in Cartesian space such that a rotation, £, about a will rotate 
frame 1 to frame 2. It is assumed that the quaternion has unit length such that 


WQaxis— angle || -L • 

In order to formulate an error signal based on quaternions, an error angle vec- 
tor based on the deviation between the command and the system state vector is 
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computed. Suppose the quaternion command is represented by q c and the system 
quaternion state is given by q. Then the quaternion representing the error between 
these two quaternion vectors is given by [63] 

QiQa. + <?2 q C2 + qzq C3 + <m C4 
- qiq C2 - <m C3 + m Ci 
qsq Cl + q 4 q C2 - qiq C3 - <?2g C4 
q±q ci - qsq C2 + q2q C3 - < m C4 

Applying a quaternion-to-Euler angle transformation for a z-axis, y-axis, and x-axis 
sequence, one can compute the this error vector in terms of Euler angles: 

(j) e = arctan 

9 e = arcsin 
ip e = arctan 

If the error between q c and q is small, the sign of q e \ can be taken as positive. In this 
case, one can apply a small angle approximation. For quaternions, this implies that 


2 (?e3?e4 + 

.2 (<?el<?el + g e 4<?e4) ~ 1 

[2 (<?el<7e3 - <7e2<? e 4 )] 

2 (? e 2?e3 + 9el9e4) 

2 (?el9el + qe2qe2) ~ 1 


q e = q c 1 ® q = 


q e i ~ l 
llfell < i 

(9.3) 

II&3II < 1 

1 1 <?e4 1 1 < 1 

This can be seen be simply examining the axis-angle formulation of a quaternion. 
With this assumption, one can approximate the error angles as 


(j) e « 2 q el q e2 

0 e ~ 2g e ig 6 3 ( 9 - 4 ) 

-0 e ~ 2g e ig e4 

This error is fed into the nominal control law along with the components of the body 
angular rate vector, u>, to compute the control signals in the roll, pitch, and yaw axes. 
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9.2.3 ifoo-NMA Architecture Formulation 


This section formulates a low-order decoupled H ^- NMA architecture for the CLV. 
The dynamics for the adaptive law design can be derived by considering the general 
attitude control problem. In this problem, the plant dynamics can be approximately 
captured as 

9 = T(q,u) 
x = cr(q, x, to, <5) 


to = f(q, x,to) + BS 

where q is the quaternion representing the inertial attitude, to is the angular velocity, 
x are the position and velocity dynamics, B is a diagonal control effectiveness matrix, 
5 is the control action, and 


9 = T(g,o;) = - 


-92 

-93 

-94 

9i 

-94 

93 

94 

9i 

-92 

-93 

92 

9i 


to = -fLcu 


(9.5) 


The control effectiveness matrix is assumed diagonal due to the restriction that the 
control law be decoupled. If one wanted to couple the adaptive control design, this 
restriction could be removed. 

The quaternion guidance command to the attitude control system is open loop [46] 
and slowly varying. Since, the quaternion is slowly varying, it is assumed that 


q c ~ 0 

This assumption allows one to formulate the CLV attitude control problem as a 
stabilization problem. To this end, the attitude dynamics are reexpressed relative to 
the slowly varying attitude command, q c . The quaternion attitude error was expressed 
in equation (9.2) as 

Qe = q^ 1 ® q 
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Since q c ~ 0, the time derivative of q e is given by 


q e = q c 1 ®q 

= <g> figca 

The quaternion attitude can be expressed as 

q = q c ®q e 


Hence, 



< 8 > 


Using quaternion algebra, 


Q Cl q ei q C 2 q e2 q C 3 q e3 Qc^qe 4 


q c ®q e = 


q Cl q e2 + q C 2 q ei + q C 3 q e4 ~ q Ci q e3 
q Cl qe 3 - q C 2 q e4 + q C 3 q ei + q C 4 q e2 


q ci q ei + q C 2 q e3 - q C 3 q e2 + q C 4 q ei 


where q Ci is the i th element of q c and q ei is the i th element of q e . From equation (9.5), 

^q c ®qe 



- (q c <8> qe ) 2 

(q c ®qe) 1 

(q c <8> g e ) 4 

- (q c <8> q e ) 3 


- (q c <8> q e ) 3 

- (q c <8> g e ) 4 

(g c ®<? e ) 1 
(9c <8) g e ) 2 


- (q c <8> g e ) 4 

(9c ® 9e) 3 

- (9c <8> q e ) 2 
(q c ®q e ) 1 


where (g c <g) g e ) i is the element of q c <g) g e . Since, g c is a unit quaternion, q c 1 is 
given by 


9d 



-9c 2 

-9c 3 

-9c 4 
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From this, after some algebra, the q e dynamics in equation (9.6) can be expressed as 

Qe2 Qe3 Qe/l 
Qel g e4 Qe% 

QeA Qel Qe 2 

~Qe 3 Qe 2 Qel 

Assuming that ca remains small and applying the small angle properties of q e in 
equation (9.3), the q e dynamics can be approximated as 

. 1 o 

ri 

q *~2 

u 

From the definition of the error angles in equation (9.4), 

0e 

6 e (9.7) 

Let the state vector for the 7 /qo-NMA architecture in figure 9.5 be defined as 

6 e 

(9.8) 

tpe 

LO 

Using this definition, the system dynamics can be approximated as 

LO 

£= (9.9) 

f c (q,x,u) + B5 

where the remaining vehicle dynamics (position and velocity) are unmodified in the 
form as previously defined as 

x = a(q, x, u>, 5) 
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With q c regarded as nearly constant, the error dynamics in equation (9.9) can be 
rewritten as 

L 0 

+ B5 (9.10) 

f c (q e ,x,u) 

where / c (-) is the equivalent form of /(•) as a function of q e instead of q and 

0 

B = 

B 

The total augmented control effort is defined as 

8{t) = 5 n (t) - 8 ad (t) (9.11) 

where 6 n (t) is the nominal control output and 5 a dit) is the augmented adaptive signal. 
It is assumed that the linear control law was designed to achieve a second order 
response in each control channel. Based on the definition in equation (9.11), it is 
assumed that the application of the nominal control law, has the effect of creating 
the following approximate form for the error dynamics in equation (9.10) 

t(t) = A m e(t) - BS a d(t) + BA(e(t)) (9-12) 

where A(e(t)) is the modeling error that exists between the desired dynamics and the 
actual dynamics and 

0 / 

A 

-k p -k d 

k p 0 0 

kp = 0 0 and K d 

0 0 k r 
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ki and b t are chosen to match the damping and stiffness associated with the desired 
second order response of the gain scheduled CLV control system design. These as- 
sumed dynamics match the form of the dynamics in Chapter 4 in equation (4.7) with 
B m = 0. Hence, one could use these dynamics to formulate a coupled H 0 0 -NMA 
architecture based on Chapter 4 for the CLV. 

In order to maintain a decoupled adaptive design, it is assumed that A (e c (t)) is a 
diagonal uncertainly of the form 


Ai (</> e ,p) 0 0 

A(e(t))= 0 A . 2 {0 e ,q) 0 

0 0 A 3 (i/) e ,r) 

This allows the dynamics of each control channel to be expressed as 


(9.13) 


Xi(t) = A m .Xi(t ) - BiS adi (t) + BiA^Xi) (9.14) 

where i represents the i th control channel, i = 1 corresponds to the roll channel, i = 2 
corresponds to the pitch channel, i = 3 corresponds to the yaw channel, X\ = [<j) e p ] T , 
x 2 = [9 e q ] T , x% = ['ipe r] T , Bj = [0 B(i,i)} T , and A mi is defined as 


A 


mi 


0 

-ki 



(9.15) 


Only the pitch and yaw channel are augmented with an adaptive control law. The 
state emulator for these adaptive laws can be thought of as an error state emulator 
because the state vector is based on the attitude command error, q e , and the angular 
velocity, c o. Note that there is no reference command in the dynamics in equation 
(9.14). This implies that, unless there is an initial attitude error and/or an angular 
rate, the reference model does not need to be implemented. In this case, the H^- 
NMA architecture simplifies. Comparing with figure 4.1, the architecture in each 
channel reduces as shown in figure 9.5. In this case, from the error definitions in 
Chapter 4, the state emulator error for the i th channel is e^(t) = Xi(t ) — Xi(t) and 
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the state emulator tracking error for the i th channel is e,(t) = —Xi(t). In this figure, 
K x represents the nominal control system for the i th channel which feeds back on the 
same error state, Xi(t). 



Figure 9.5: H ^- NMA architecture simulation diagram for the CLV. 


9.3 Simulation Results 

The adaptive system and the nominal control system is implemented in discrete time 
with an update rate of 50 Hz. All system latency and actuator limitations are mod- 
elled. The actuator command from the control system is in terms of a roll RCS 
command, a pitch angle command to the gimbal, and a yaw angle command to the 
gimbal. However, the actual gimbal control signal is a tilt and rock command. The 
tilt and rock command is computed based on the following equation. 


^ tilt 


cos(4>tr ) 

sin((f>TR ) 


fipitch 

^ rock 


- sin((pTR ) 

cos(4>tr ) 


3yaw 


where 4>tr is rotation angle of the tilt-rock gimbal, Utut is the tilt angle command, 
Urock is the rock angle command, 8 rntc h is the control system gimbal pitch command, 
and 5 yaw is the control system gimbal yaw command. 
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The iJoo-NMA design was meant to be minimum complexity. Therefore the vector 


/3(xi) was chosen as 

a 

P(Xi) = 

Xi 

where a is a user defined constant and Xi is the corresponding state vector of the i th 
channel in equation (9.14). It was assumed that the desired response is equivalent to 
choosing a damping ratio of Q — 0.707 and an undamped natural frequency of u n = 5 
rad/sec in both the pitch and yaw control channel. This implies that the stiffness 
and damping parameters of A mi should be defined by 

ki = u 2 n and 6* = 2 (co n 

The adaptive gain was simply set to T = 30. The linear H 0 0 design was augmented 
with an integrator. The augmented state vector used for each design was 

Xinti 

X = 

Xi 

where x inti is the integrator state. Using this augmented state vector, the design 
strategy in section 4.4 was followed. In each design, the design matrices Q 5 and R 5 
were chosen as 

Q* = I and R? = 0.025 

This selection of weighting matrices suggests that the tilt-rock actuator should have 
at least a 3 Hz bandwidth. 

The simulation results focus on the mean international space station mission and 
examine the degree to which the nominal control design can be improved by adding 
a decoupled state feedback adaptive element described in the previous sections. To 
measure this degree of improvement, tracking performance is compared between the 
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baseline control performance and the H ^- NMA architecture performance for the 
Worst- on- Worst( WoW) Monte-Carlo dispersion cases. The WoW cases capture the 
combination of the worst possible uncertainties^. e. dispersions) occurring simulta- 
neously. Simulations show that the adaptive control law always improves the perfor- 
mance of all of the stable WoW cases. Results show that system tracking errors and 
measures of structural stress are reduced. In many examples of adaptive control law 
implementation, the control laws exhibit increased control effort. This could be detri- 
mental to the performance of a launch vehicle. However, in all of the WoW cases, the 
total control activity between the adaptive cases and the linear control cases remain 
approximately the same (as measured by duty cycle and duty cycle rate). Moreover, 
in the case of rock duty cycle rate, control activity is generally reduced. For the 
WoW cases that are unstable with only the linear control law, the adaptive control 
algorithm is able to maintain acceptable tracking performance in almost all of the 
cases while control effort remains consistent with other dispersion cases. 

In comparing performance of the adaptive control system algorithm with the base- 
line flight control system, parameters that characterize attitude tracking, use of ef- 
fectors, loads, and errors at the time of first stage separation were captured from the 
simulations and evaluated. Performance emphasis is on loads and attitude control. 
These metrics are listed below. 

• Total Q-Alpha - Square root of the sum of Q-Alpha (aerodynamic pressure 
multiplied by angle of attack) and Q-Beta (aerodynamic pressure multiplied by 
side slip angle). Reduction of this is desirable in all cases since it represents 
aerodynamic loading on the vehicle, therefore it is heavily weighted. 

• Attitude errors - roll, pitch, and yaw errors (command minus sensed) before 
filtering. This is coupled with total Q-Alpha during high Q since guidance is 
essentially commanding zero aerodynamic angles in this region of flight. Due 
to thrust vector dispersions, pitch and yaw couple into roll and affects RCS 
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propellant consumed. 


• Total nozzle gimbal angle - The maximum value should be maintained below 
the specified value due to hardware capability. This metric is heavily weighted 
since exceeding capability could result in loss of the vehicle. 


Gimbal duty cycle - area under the total nozzle angle curve as defined by 


Duty Cycle (t) = / (Nozzle Gimbal Angle) dt 


Jo 

This is not heavily weighted, recognizing the fact that utilizing the effectors 
more aggressively may be necessary to achieve better results. 


• Gimbal rate duty cycle - number of sign changes in both the rock and tilt 
actuator rates. This metric is needed to evaluate actuator chatter. It is not 
heavily weighted. 


• Body rates (truth, not measured) at separation. Roll rate is not as important 
as pitch and yaw. Large rates can cause interference between first stage and 
interstage hardware. Pitch and yaw are weighted more heavily than roll. 


To capture these factors, a scoring metric was developed. Success was judged 
based upon increasing the maximum values of the performance metrics. Simulations 
were made with the baseline flight control system and the adaptive control system 
using the same set of dispersions. For each simulation, a scaling, S(i), for each metric 
is computed. This value is given by 


S(i) 


1 + 


B(i) - M(i ) 

W) 


( 9 . 16 ) 


This expression is based on a percentage difference in the bounding values, B(i), 
and the maximum value of each metric during a simulation run, M(i). Below are 
bounding values and the weights for the performance metrics used to score each of 
the simulations (values were suggested by NASA). 
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The score for each metric is 


Score(i) = W{i)S(i) (9-17) 

and the total score for each simulation is 

Total Score = Score (i) (9.18) 

i 

The metric score has improved if the total score increases. 

The i^oo-NMA adaptive control architecture improves the metric score for all of 
the Worst-on- Worst cases in which the baseline control law is stable. Figures 9.6 
and 9.7 show the raw metric scores for the LGo-NMA architecture and the baseline 
control law. There are 100 WoW cases. Case numbers without data represent a 
case where the baseline control law is unstable (no comparison can be made). A 
case is considered unstable if the roll error exceeds 10°, the pitch error exceeds 5°, 
or the yaw error exceeds 5°. For the baseline linear control law, WoW cases 7, 10, 
46, 47, and 82 are unstable. For the adaptive control law, only cases 46 and 82 are 
unstable. Figure 9.8 shows the percentage improvement in the performance metric 


Table 9.2: CLV Metric Bounding Values 


Metric, M (i) 

Bounding Value, B(i) 

Total Q-alpha 

Attitude errors (absolute values) 
Total nozzle gimbal angle 
Gimbal duty cycle 
Gimbal rate duty cycle 
Body rates at separation (absolute values) 

5500 2^ 

Roll-20 0 , Pitch-3°° and Yaw-3° 
4° 

193 -2- 

sec 

Rock-30 cycles, Tilt-30 cycles 
Roll-60 — Pitch-5 — Yaw-5 — 

see 7 see 7 see 


Table 9.3: CLV Metric Weights 


Metric, M (i) 

Weight, W{%) 

Total Q-alpha 

5 

Attitude errors (absolute values) 

Roll-2, Pitch-3, and Yaw-3 

Total nozzle gimbal angle 

5 

Gimbal duty cycle 

2 

Gimbal rate duty cycle 

Rock-2, Tilt-2 

Body rates at separation (absolute values) 

Roll-2, Pitch-5, Yaw-5 
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for the adaptive control system compared to the corresponding baseline control case. 
Figure 9.9 shows the percentage improvement in the performance metric when the 
first 20 seconds is excluded from the metric calculation. 

Next the time history of each WoW dispersion is compared. In all of the following 
hugres, the unstable cases for both the adaptive and baseline control law are not 
shown for clarity. Figures 9.10 - 9.12 show the attitude error for each stable baseline 
and adaptive control case. In each case, roll error between each baseline and adaptive 
case stays about the same but this isn’t surprising because the roll control law is 
not augmented with an adaptive element. The peaks in the pitch and yaw error is 
generally reduced using the adaptive controller. Figures 9.13 and 9.14 show the angle 
of attack and sideslip respectively. These plots show that peaks in both quantities 
are generally reduced. Figures 9.15 - 9.17 show plots of variables related to structural 
stress. Once again, in each case, peaks are reduced. Figures 9.18 and 9.19 show the 
rock and tilt command respectively. In these figures, the total control effort remains 
approximately the same. This is good because in many applications of adaptive 
control, control activity tends to increase relative to the linear baseline design. Figures 
9.20 - 9.23 show the duty cycle and duty cycle rate for the tilt and rock actuators 
during each simulation run. While the measures of duty cycle remain approximately 
the same, it is notable that the tilt and rock duty cycle rate decreases as a general 
trend. 

9.4 Conclusions 

An T/qo-NMA architecture for the Crew Launch Vehicle was developed in a state 
feedback setting. The minimal complexity adaptive law was shown to improve base 
line performance relative to a performance metric based on Crew Launch Vehicle 
design requirements for all most all of the Worst-on- Worst dispersion cases. The 
adaptive law was able to maintain stability for some dispersions that are unstable 
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with the nominal control law. Due to the nature of the H ^- NMA architecture, the 
augmented adaptive control signal has low bandwidth which is a great benefit for a 
manned launch vehicle. 
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WoW Case Number 


Figure 9.6: Raw metric score values for the stable baseline control WoW cases. 
Black is the adaptive control law and grey is the baseline control law. 
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Figure 9.7: Raw metric score values after 20 seconds. Black is the adaptive control 
law and grey is the baseline control law. 
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Figure 9.8: Summary of metric score improvements for stable WoW cases. 
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Figure 9.9: Summary of metric score improvements for stable WoW cases after the 
first 20 seconds. 
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Figure 9.10: CLV roll error comparison. Black lines represent the baseline responses 
and red lines represent the H ^- NMA responses. 
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Figure 9.11: CLV pitch error comparison. Black lines represent the baseline re- 
sponses and red lines represent the i/oo-NMA responses. 
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Figure 9.12: CLV yaw error comparison. Black lines represent the baseline responses 
and red lines represent the H ^- NMA responses. 



Figure 9.13: CLV angle of attack comparison. Black lines represent the baseline 
responses and red lines represent the L/qo-NMA responses. 
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Figure 9.14: CLV sideslip comparison. Black lines represent the baseline responses 
and red lines represent the H^-NMA responses. 
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Figure 9.15: CLV Q — a comparison. Black lines represent the baseline responses 
and red lines represent the H ^- NMA responses. 
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Figure 9.16: CLV Q — /3 comparison. Black lines represent the baseline responses 
and red lines represent the H^-NMA responses. 
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Figure 9.17: CLV Q — a total comparison. Black lines represent the baseline 

responses and red lines represent the L/oo-NMA responses. 
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Figure 9.18: CLV rock command comparison. Black lines represent the baseline 
responses and red lines represent the H^-NMA responses. 
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Figure 9.19: CLV tilt command comparison. Black lines represent the baseline 
responses and red lines represent the L/qo-NMA responses. 
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Figure 9.20: CLV rock duty cycle comparison. Black lines represent the baseline 
responses and red lines represent the H^-NMA responses. 
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Figure 9.21: CLV tilt duty cycle comparison. Black lines represent the baseline 
responses and red lines represent the L/qo-NMA responses. 
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Figure 9.22: CLV rock duty cycle rate comparison. Black lines represent the baseline 
responses and red lines represent the H 0 C -NMA responses. 



Figure 9.23: CLV tilt duty cycle rate comparison. Black lines represent the baseline 
responses and red lines represent the H^-NMA responses. 
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